Meson bound states in multiflavour massive 



Schwinger model* 



M. SADZIKOWSKI t 

Institute of Nuclear Physics, 
Radzikowskiego 152, PL-31 342 Krakow, Poland 

and 

P. Wegrzyn* 
Institute of Physics, Jagellonian University, 
ul. Reymonta 4, PL-30 059 Krakow, Poland 



Abstract 

The problem of meson bound states with Nf massive fermions in 
two dimensional quantum electrodynamics is discussed. We speculate 
about the spectrum of the lightest particles by means of the effective 
semiclassical description. In particular, the systems of fundamental 
fermions with SU(2) and SU(3) flavour symmetries broken by massive 
terms are investigated. 



*Work supported by the KBN grants no. 2P 30207607 and 2P 30204905. 
^Fellow of the Polish Foundation (FNP) scholarship for the year 1996 
* E-mail wegrzyn@ztc386a.if.uj.edu.pl 



1 



1 Introduction 



The literature on the Schwinger model |T| and its generalizations should 
refer to at least a hundred of interesting and important papers. This the- 
oretical model is usually used to demonstrate some important phenomena 
known from more realistic models. It is easy to see and investigate here 
screening and confinement, U(l) problem, chiral symmetry breaking and 
topological vacua, formation of bound states and finite temperature effects. 
On the other hand, known non-perturbative techniques tested here, like lat- 
tice calculations or sum rules, reveal many pitfalls, what makes their use 
doubtful. Surprisingly, the two-dimensional Schwinger model shares many 
phenomenological features in common with four-dimensional QCD. While 
one can state that the problem of massless charged fermions in two dimen- 
sions has been elaborated in all details the situation when fermion 
masses are finite [12]-pl|] still requires further investigations. In particular, 
it would be interesting to understand in the theory with the confinement 
the mechanism of formation of the lightest physical particles from the fun- 
damental fermions. This subject is discussed in this paper. 

Our paper is organized as follows: in the Section 2, we introduce briefly 
the multiflavour massive Schwinger model in its bosonized version. Classi- 
cal equations of motion and ground states are described. In the Section 3, 
we discuss the lightest meson bound states in the case when fundamental 
fermions ('quarks') are heavy (or in other words the coupling constant is 
weak). To make our paper self-contained, in the first part of this section we 
remind the most important points of the semiclassical quantization proce- 
dure applied to the 'particle-like' solutions of the sine-Gordon theory. This 
knowledge is used in subsequent developments. In the Section 4, the analy- 
sis of the lightest meson states is extended to the case of light quarks (the 
strong coupling regime). We discuss the model for SU(2) and SU(3) flavour 
groups separately, both in unbroken and broken cases. In the Appendix, we 
advertise the method useful to derive approximately particle-like solutions 
to the classical field equations. This method is explained on the simplest 
example of the sine-Gordon theory, but its advantage is that it can be used 
to more involved problems. 
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2 Lagrangian 



We will consider Nf fundamental charged and massive Dirac fermions in 
(l+l)-dimensional Minkowski spacetime, with the following Lagrangian den- 
sity 

1 N f 

C = —F^ v + frii-fdr - vfAp - m a )ip a , (1) 

a=l 

where 

= d»A v - dyAp, (2) 

e is the coupling constant (dimension of mass). We allow for masses of 
fermions m a to be different for each flavour a. In the exactly solvable mass- 
less case m a = 0, the classical system possesses a global symmetry group 
U{N) L x U(N) R = U{l) v x U{l) A x SU{N) L x SU(N) R . At the quantum 
level, the axial symmetry U(1)a is broken down by the anomaly. In the 
massive case, the flavour symmetry SU(N)l x SU(N)r is broken explicitly 
to U(l) N f~ 1 , leaving only (Nf — 1) conserved flavour numbers. We are going 
to make use of the standard (Abelian) bosonization rules, 

Nmairrr] = -^e^d u $ a , (3) 



N ma bP a ^ a ] = -cm a N ma [cos V^$ a ] , (4) 

where N ma denotes normal-ordering with respect to the fermion mass m a , 
$ a is the family of canonical pseudoscalar fields. The appearance of the 
constant c = e 7 /27r (7 is the Euler constant) is due to the use of the specific 
renormalization scheme. Obviously, physical quantities do not depend on 
this choice. Then, one can derive the bosonized version of ([l]): 

Nf 

+ ^2 cm l N m a [cos V4Tr$ a ] + const. (5) 

a=l 

The different vacua are labelled by the angle parametr 6 and the relevant 
constant in (J5[) should be adjusted to ensure zero vacuum energy. After 
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integrating out the electric field Fqi, we arrive at the effective Lagrangian 
describing the system of interacting pseudoscalar fields defined by 

a=l L 



2 / JV / 



N f „ \ 2 N 



v -ff = 2^ ( E * a + ^ J - E ^ [ cos ^* a l + const • ( 6 ) 



The topological charges of pseudoscalar fields <£ a are related to flavour quan- 
tum numbers of fundamental fermions through the relations: 



Q a = 

V7T 



+oo 



Classical Euler-Lagrange field equations can be derived, 

o2 



(7) 



d^ a - d 2 x § a + — \ Y & + -j= + cv^rm* sin (V^$> a ) = . (8) 

7T \ V47T / V ' 



v 6=l 

Classical vacua are easy to determine from the requirements 
dV eff /d$ a = 0, namely: 

I A Nf 

$ a = -= arcsin — ^ + V^n" , n a G Z , V n a = , (9) 



; — arcsin — » + Jnn a , ra a G Z , > 
where the constant ^4 is subject to the equation 



A 2-7T 

V arcsin +6* + c(—) 2 A = 0. (10) 
, mi e 

a=l a 



All finite-energy (localized) solutions should approach asymptotically the 
vacuum: 

^4tt 



<3? — ► — = arcsin — ^ + y/Tml , (11) 



d^ a , (12) 

E*4 = 0. (13) 
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There are no restrictions on our considerations if we set all n a _ equal zero. 
Thus, we have n°\ = Q a , and Eq.(|l3|) means that all finite-energy solutions 
are chargeless (charge screening). 

For 9 = (no CP breaking) we find 

N f 

$ a vac = ^n\ n a eZ, ]>> a = > ( 14 ) 

a=l 

and we fix the constant in (JgJ) to obtain the potential 

N f 2 /N f \ 2 

Veff = E Cm l N m a [1 " COS V4^$°] + ^ £ > ( 15 ) 

a=l \a=l I 

For 9 7^ 0, the equation ( |Toj ) has still a solution corresponding to a 
classical vacuum. However, the case of 9 = ±7r is special, and the degeneracy 
of the vacuum structure appears. 



3 Heavy quarks. 

At first, we consider the case when the Lagrangian (bare) masses of funda- 
mental fermions are much larger than the scale of electromagnetic interac- 
tions e. 

Let us divide the Lagrangian @ into two parts, 

Nf Nf ^ 

4 = E2 (w2 + ^* jcos ^ ri ' (16) 

a=l a=l 

P 2 f N f ft V 

^JH-SF^ — ^) ■ ("> 

The first part defines the system of Nf sine-Gordon fields. The weak in- 
teractions between them ( |l7|) are important when fields are close to their 
vacuum values. 

Since the theory is superrenormalizable, in order to subtract all ultravi- 
olet divergencies it is enough to replace the unordered functions of fields in 
the Lagrangian by their normal-ordered counterparts. Moreover, we do not 
need to take much care of the renormalization of quadratic forms of fields 
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(kinetic terms and Coulomb interactions), because it effects only in the ad- 
dition of an infinite constant to the Lagrangian. Thus, everything we want 
to know about the ultraviolet renormalization is contained in the prominent 
formula derived by Coleman [pffl ]: 

cos (/?$)=( -J N m [cos (/?$)] , (18) 

where A is the ultraviolet cutoff. The above formula allows us to compare 
two different renormalization scales, 



N m 



[cos (£$)] = (£) ^ [cos (/?$)] . (19) 

\m J 



Thus, the effect of finite renormalization is a multiplication by the power of 
mass ratios. For free fermions (|T^), the anomalous dimension (3 2 /Att is equal 
one, thus the operator ( |l8| ) acts just like a free fermion mass operator. 

If we want to have the standard form of sine-Gordon Lagrangian, we 
need to renormal-order CP e ^ (|l6|), 



-eff 

o=l 



{d^ a f + ^ cos ((3$ a ) 



(20) 



2 V /x , p 

where (3 = y/4iv and the renormalized mass is: 

M a = (3 2 cm a = 2e^m a . (21) 

As the system of classical field equations derived from C®ff + ^jff ( 1Q > 
17) admits localized (particle-like) solutions, the semiclassical quantization 
can be used to find, at least the lowest, bound states. 

There exist two different types of particle-like solutions for the sine- 
Gordon equation. The first one refers to solutions which are time-independent 
in their rest frames (static, solitonic), called usually solitons (Q a = +1) and 
antisolitons (Q a = — 1), 

$ a = Q«l arctan ( e -M»(*-*»)^ . (22) 

The energy density for such a solution is localized around some point x a . The 
classical masses are M% = 8M a /(3 2 . As the quantum (finite) corrections 
do not introduce any new mass scale, we can denote the quantum mass of 
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solitary solitons as M. a qu = 8M a /(3' 2 , and quantum effects are reduced to 
the change of the coupling constant. The semiclassical WKB quantization 
around solitons or antisolitons gives the following effective coupling constant 



,/2 _ P P=V^ 



r = ^r^ Svr. (23) 



1 - 2- 

1 8tt 



For f3 = \/47r, the same result has been obtained by the calculation based on 
Feynman integrals p3| and it is very close to the numerical value obtained 



with variational methods ^6[. Presumably, for (3 = v4tt the result of WKB 
approximation can be expected to be exact. 

The second type of particle-like solutions for the sine-Gordon equation 
corresponds to periodic in time solutions, called breather modes: 

4 / 7 smut \ , nd . 

$ a = - arctan J- , (24) 

p \u) cosh 7a; / 



where T = 2n/to is a period and 7 = yj (M a ) 2 — lo 2 . These solutions are 
localized and 'topologically chargeless' Q a = 0. Their quantization via WKB 



methods [23, 24 1 is analagous to the quantization of the Bohr orbits of the 



hydrogen atom. The Bohr-Sommerfeld quantization condition is here: 

S + ET = 2irn , (25) 

where n is positive integer, S is the action per one period, and E is the 
energy. Again, it can be checked that all quantum effects can be acknowl- 
edged in the change of coupling constant according to ( p3| ) . The masses of 
quantum states produced on the basis of classical breather modes are: 

l6M a n(3' 2 . . 

• M - = ^ sin ^"' (26) 



where n = 1, 2, 3, ... < 8tt/P' 2 . For j3 = \[Ar there are no such states. 

After this brief description of the spectrum of the quantized sine-Gordon 
system, we can pass to the discussion of the quantum states for the Schwinger 
model (16,17). The main result of interactions ( |T7| ) is to impose the zero 



total charge condition (|l3l). The lightest bound states is composed of a pair 
of a soliton and an antisoliton. Its mass (up to e/m a corrections) is given 
by 



2p 7 

M=M a qu + M b qu = —(m a + m b ) . (27) 
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Let us specify 6 = and M a = M b = M, and Q a = -Q b = 1. Then, the 
corresponding pair of wave packets (|^) is an exact solution of equations 
of motion provided that x a = Xf,, i.e. the soliton and the antisoliton are 
localized at the same point. It can be immediately checked that the energy 
of Coulomb interactions ( |T7| ) is zero, so that both particle- like waves do 
not interact together. The state is composed of a free soliton and a free 
antisoliton, living together around the same point x a = Xf,. But if we try 
to separate them for a small distance Ax = x a — x^, they start to interact 
and the energy increase due to Coulomb interactions (|l7]) can be easily 
calculated, 

AE = ^f(Ax) 2 . (28) 

We have observed that the asymptotic freedom and the confinement in the 
Schwinger model can be already seen at the classical level, provided that the 
bosonization has been performed. In this way, the bosonization provides us 
with some kind of the dual description of the model. 

The meson state described above were composed of a pair of wave packets 
( f22"l) corresponding to two particle-like solutions of different flavours. There 
exist also states where 'quark' and 'antiquark' wave packets are of the same 
flavour. However, these solutions are of different type, being breather modes 
for the system defined by the equations (||). Their masses (close to 2M a ) 
and their shapes can be found using the approximate method described in 
the Appendix. 



4 Light quarks 

In the strong coupling limit, where e 3> m a it is much more convenient to 
perform appropriate change of field variables: 

x a = 0^ + -=£=<5? (29) 
^/47riV j 

using the orthogonal matrix O defined below: 



O 



^N f (N f -l 



= (1,1,... 1,-^ + 1) 
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o Nf 



y/(Nf-l)(Nf-2) 

-L(l,-1,0,...,0) 



(1,1,... l,-iV> + 2,0) 



(30) 



Using new fields the Lagrangian density (||) takes the form: 



£eff = -a Y.^f-y&f+Y, cmlN ma [co S V^(0 T xT - w ]+const 



a=l 



0=1 



(31) 

where /i 2 = Nje 2 /it. The topological charges are defined now with respect 
to x a fields, and they can be interpreted as the electromagnetic charge 



Q 



x 



+oo 



(32) 



and SU(Nf) hypercharges for other fields. For instance, in the case of two 



flavours N 



/ 



X 



(33) 



is a third component of isospin. For Nf = 3 we have two conserved hyper- 
charges, 

2 1 



Y 



i . 



X 



X 



+oo 
— OO ; 
+ 0O 
— OO 



(34) 



(35) 



Following Coleman fl~3[| , we can decouple the only heavy field x 1 . It is 
straightforward to perform necessary renormal-ordering, 



2V„ 




4tt ! v 



and remove x 1 from (31) 



flight 
^eff 



m. 




4tt. i ^ . 



(36) 



a=2 



9 



2 ( V> \ N f 



+ cm 2 a I 



a=l 



9 



N 



f 



N f 6=2 



+ const, (37) 



The simplest case Nf = 2 was already elaborated by Coleman [13|]. He 
noticed that the Lagrangian density (37) is equivalent (for Nf = 2) to the 
sine-Gordon theory, 



4fi = \M? + ^ M " N m [cos (v^x 2 ) 



where 



M 



e' a 1 ! 2 \jm\ + mf. + 2m\m2 cos 



2/3 



(38) 
(39) 



The quantum states of the sine-Gordon theory were desribed briefly in the 
previous section. Here we have (3 = V2n, what implies that there are two 
breathers (|2~4]) with masses Ai and (A4 is here the quantum mass 

of the soliton). Coleman identified three solutions: the soliton [Q = 0, 
^3 = +1), the antisoliton {Q = 0, I3 = —1) and the lighter breather {Q = 0, 
I3 = 0) as the components of the isotriplet. These are the lightest physical 
states. Coleman noticed that they form the degenerate isotriplet even if the 
SU{2) flavour symmetry is apparently broken (i.e. mxjmi is very large or 
very small). However, it is not true that the SU(2) symmetry is restored 
here exactly. The differences in the masses of the isotriplet states arise as we 
take into account the corrections of the order m a /u. Note that the origin of 
J3 = state (the analogue of ir°) is different that the origin of J3 = ±1 states 
(the analogues of ir^). Another interesting thing is that the dependence on 
the vacuum parameter 9 is only via the sine-Gordon mass value. As far as 
the heavier breather is concerned (that of mass A4^/3), Coleman wrongly 
|L9f identified that as the isosinglet state (the analogue of r\ particle). In fact, 
the isosinglet state should be matched with the non-trivial configuration for 
Xi, so that its mass is of the order u. Then, the heavier breather mode 
should be rather interpreted as a bound state of pions (or ir° excitation). 

Now, let us consider Nf = 3 case. We restrict ourselves to the case when 
the vacuum angle 9 is zero (C, P symmetries are not broken). Our effective 
Lagrangian contains now two light fields, 



C 



light 
N f =3 



1 



2\2 



v, 



eff . 



eff 



c» 1/3 m[ /3 N mi 



[2tt 

cos (V2^x 3 + y yX 2 ) 



+ 
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ml /3 N m2 



i- cos (\/2^x 3 -y T x 2 ) 



5/3 

*3 1V "13 



1 — cos \ —x 



(40) 

It gives the following field equations of motion (a = 2,3), 

r)V cc 

dh a ~d 2 x X a + -^f = 0- (41) 

At first, we discuss the case when all fermion masses are equal, namely 
mi = m<i = 1TL3 = m, so that SU(3) flavour symmetry remains unbroken. 
We list several exact classical solutions, which correspond to the lightest 
bound states. At the begining, note that if the field x 2 takes its vacuum 
value, say x 2 = 0j the equations (H) reduce to the sine-Gordon equation 
for the field x = X 3 > 

M 2 

d 2 x-d 2 xX + ^=sin(V2^x) = , (42) 

V Z7T 



where 

M= (2e 7 /x 1/3 m 5/3 ) V2 . (43) 

Therefore, a soliton, an antisoliton and the lighter one of breathers represent 
here three solutions with equal masses, and quantum numbers Y = and 
^3 = +1,-1,0 respectively ('pions'). Because of the exact SU(3) flavour 
symmetry, there are more solutions of the same mass as pions. In order to 
obtain them, we need to put x 2 /V3±x 3 = an d observe that x = xVv^T 
X 3 satisfies the equation ( |42"| ) (the sign corresponds to two alternatives). 
Solitons an antisolitons constructed in these two cases describe four bound 
states with quantum numbers Y = ±1 and ^3 = ±1/2 ('kaons'). On the 
other hand, periodic solutions (breathers) for both alternative cases can 
be identified, so that the lighter breather gives us the eighth state of the 
pion mass ('rfe-particle'). In this way, we have completed the whole meson 
octet, being the family of the lightest physical particles. Let us remark that 
this symmetric case allows to construct further exact solutions. Heavier 
breathers give next three particles (mass/octet mass= VS,Y = 0,±1, / = 
0). Another example we obtain if we assume X 3 = and construct the 
soliton solution for the corresponding equation for x 2 (mass/octet mass = 
3/2 + \/31og (5 + Y = ±2 and / = 0). 

We now turn to the discussion of the case when m = mi = 1112 and 
[i S> 77T-3 ^> m, i.e. the 'strange' quark is much heavier but still below the 
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scale of interactions. The states of the lowest isotriplet satisfy the equa- 



tion (42) with the same assignment of the solutions as before. The lowest 
isodoublet states can be constructed in the similar way as in the case of 
equal masses. The only difference is that the combinations of the fields 
X 2 /\/3 =t x 3 do not take vacuum values, but they remain still small with 
respect to x = X 2 /v / 3T X 3 - Solving the static equations of motion (|4l"|) 
(using methods given in the Appendix) for these variables we found that 
the mass ratio between isotriplet states (pions) and isodublet states (kaons) 
is approximately 7r 2 y / r/4v / 3 where r = (777,3/m) 5 / 3 . 

Finally, we discuss the case when m = m\ = ni2 and 771,3 ^> fJ- ^> m, i.e. 
one of the flavours refers to heavy quarks. In this case we rotate both of 
the light fields & a (a = 1,2) through the transformation (29) into the new 



field variables x a > and heavy field <3? 3 remains here as the third field variable. 
Using the new fields, our effective lagrangian (^) (with Nj = 3 and 6 = 0) 
takes the form: 

A// = lid.x 1 ) 2 + \{d, x 2 ? + \{d^f - v eff 

^// = 7(x 1 + ^j - cm 2 N m [cos {x 1 +X 2 )+ 

cos v / 27r(x 1 - X 2 )] ~ cm%M mz cos V4tt$ 3 + const. (44) 

To consider the lowest states, let us put the heavy field in its vacuum state 
$ 3 = 0. Then, we derive from the equations of motion that field x 1 1S m its 
vacuum state as well, and field x 2 satisfies the sine-Gordon equation with 
the mass M = \/4tt cm and (3 = \^2tt. The three solutions of the equal mass 
to this equation, a soliton, an antisoliton and one of the breathers represent 
the three partners of the isotriplet (as it was described after Eq.(|42"|)). The 
other light states with a nontrivial contribution from the heavy field can 
also be constructed. These are isodublets I 3 = ±1/2 with nonzero flavour 
number (7) Q 3 = ±1, built on the static solutions of the equation of motion. 
Since the flavour symmetry is broken, the mass of this state differs from that 
of the isotriplet. We found the approximate value of the mass ratio isotriplet 
mass/isodoublet mass = 7r 2 7773/(12777.). We also calculated the mass of the 
77-particle (the appropriate solution is based on the breather), and its value 
calculated up to the first order of the method described in the Appendix, 
is consistent with the value predicted by Gell-Mann-Okubo mass formula 
7T7 2 + 3m 2 = 4777^ with the accuracy of 10 per cent. 
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5 Summary 



In this paper, we have described the lightest physical meson states formed 
from the fundamental fermion fields in the Schwinger model. Of course, 
the picture is dependent of the hierarchy between the quark masses and 
the scale of interactions. For heavy quarks, we have noticed that both 
confinement and asymptotical freedom can be anticipated at the classical 
level. Physical quantities are analytical both in the coupling constant and 
in the inverses of quark masses, so that both parameters can be used to define 
perturbative expansions. To describe the physics of light quarks governed by 
the Schwinger model, we find it convenient to use a different set of bosonic 
field variables. If all light quark masses are equal (i.e. SU (N) flavour 
symmetry occurs), then we have checked for SU(2) and SU(3) examples 
that the lightest physical states belong to the mesonic SU(N) multiplets, 
excluding always the singlet state. The mass of the singlet state lies near 
the scale of interactions (i.e. the Schwinger mass). It is important to note 
here that the members of the same multiplet correspond usually to several 
different types of classical solutions of the bosonized model. In the case of 
SU(2) symmetry, the breaking of flavour symmetry can be noticed only as 
a 'hyperfine' splitting of the lightest multiplet ('hyperfine' means here being 
of the order of the inverse of the Schwinger mass). In the case of SU(3) 
symmetry (and higher groups), the situation is different. When the flavour 
symmetry is broken, we notice immediately within the lightest multiplet the 
splitting of states corresponding to different isospins. 

The Appendix 

We describe here the general approximate method to derive particle- 
like solutions for the given system of non-linear field equations. We assume 
that we are dealing with two types of particle-like solutions: static solutions 
(solitons) and periodic solutions (breathers). We present the method looking 
at the simple example of the exactly solvable sine-Gordon equation, but 
there no obstacles if put it into use in more complicated cases. 

The sine-Gordon equation reads, 

M 2 

d?<$> -d 2 x <$> + — sin (J3&) = . (45) 

First, we are attempting to find a static solution localized around some 
point xq. We specify the asymptotic conditions <&{x = — oo) = and 
$(x = +oo) = 2ir/f3. Far away from the localization point xq the field 
$ is close to its vacuum value, and the equation (|45|) can be linearized. We 
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write down the linear differential equations for x « xq and for x » x$ 
respectively. Having imposed proper asymptotic conditions, we glue both 
solutions together at the point xq. The final result yields, 

_ j * e M(x-x ) if x < XQ 

®-{ 2n _ n e -M(x- X0 ) if x > Xq ^ 

It gives the mass 7r 2 M//3 2 , being close to the exact result 8M/(3 2 . Using 
this procedure, one can calculate further corrections. 

To find periodic solutions (breather modes) we need to use some more 
tricky procedure. Any periodic solution can be expanded in the Fourier 
series, 

oo 

<3? = ^2 4>n(x) sin (ami) , (47) 

71=1 

where oj = 2ir/T. We restrict here to solutions antisymmetric in time, 
<&(— t, x) = —$>(t,x). It allows us also to write down the classical energy in 
the following way: 

I r+oo 



E = - I dx 

-oo 



t=0 



2 

The Bohr-Sommerfeld quantization condition ( |25| ) reads, 

-oo r T 



(48) 



/+oo pi 
dx / dt [d t $} 2 = 2vriV . (49) 
-oo J 



The quantum effects will come to drive the coupling constant (p3|). To 
find a classical solution being the starting point for the WKB quantization 
procedure J23|, as a first approximation we take only the first term in the 
Fourier expansion fl47p. The lowest frequency u is derived from (^) (for 
N = l), 

p+oo 

oj / dx<pi(x) = 2 , (50) 



where the Fourier coefficient <j)\ is calculated from some ordinary differential 
equation, which is dependent of the parameter oj. The energy (mass) ( j48| ) 
corresponding to the lowest breather mode in this first approximation is 
given by the following simple formula, 

E = oj. (51) 

This is nothing else but the famous de Broglie relation. Here, this approxi- 
mate formula is as well as the harmonic approximation of the breather mode 



14 



solution is good enough. If we take the exact solution, we can verify that 
E/uj = 2\/3/vr = 1.102, so that it is pretty accurate. But it is important to 
stress that the above method allows to calculate further corrections if one 
wishes. 
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